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ABSTRACT 

We present interferometric observations of three ultraluminous infrared galaxies (ULIRGs; UGC 
5101, Mrk 273, and IRAS 17208—0014) in the 3- mm wavelength range, using the Nobeyama Millimeter 
Array. Both the HCN (J=1~0) and HCO"*" (J=l-0) molecular lines were observed simultaneously. 
HCN emission was clearly detected at the nuclear positions of these ULIRGs, and HCO+ emission 
was detected at the nuclear positions of UGC 5101 and IRAS 17208-0014. The HCN to HCO+ 
brightness-temperature ratios toward the nuclei of the three ULIRGs were derived and compared 
with those of lower luminosity galaxies known to be dominated by active galactic nuclei (AGNs) or 
starbursts. In UGC 5101 and Mrk 273, where there is evidence for obscured AGNs from previous 
observations at other wavelengths, we found high HCN/HCO+ ratios (>1.8) that are in the range 
found for AGN-dominated galaxies. In IRAS 17208—0014, where the presence of a powerful obscured 
AGN has been unclear, the ratio ('^1.7) is in between the observed values for starburst- and AGN- 
dominated galaxies. The high HCN/HCO"*" brightness-temperature ratios in UGC 5101 and Mrk 273 
could be the consequence of an HCN abundance enhancement, which is expected from chemical effects 
of the central X-ray emitting AGN on the surrounding dense molecular gas. Our proposed millimeter 
interferometric method based on HCN/HCO"*" ratios may be an effective tool for unveiling elusive 
buried AGNs at the cores of ULIRGs, especially because of the negligible dust extinction at these 
wavelengths. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: ISM — radio lines: galaxies — 
galaxies: individual (UGC 5101, Mrk 273, and IRAS 17208-0014) 



1. INTRODUCTION 

Ultraluminous infrared galaxies (ULIRGs), discovered 
during the IRAS all-sky survey, radiate the bulk of their 
very large lu minosities (L > 10^^L,n) a s dust emission in 
the infrared ijSanders fc Mir abelll 1996(1 . Hence, powerful 
energy sources, such as starbursts or active galactic nu- 
clei (AGNs), must be present hidden behind dust. Since 
distant ULIRGs dominate the cosmic infrared back- 
ground emissi on, the sum of the dust-obscured activity in 
the universe l)Blain et alJll999|) . establishing the energy 
sources of ULIRGs is closely related to understanding 
the connections between obscured AGNs and starbursts. 
Distant ULIRGs are too faint to investigate in detail with 
existing facilities and so the study of nearby ULIRGs 
continues to be an important means of understanding 
the distant ULIRG population. 

It is now generally accepted that the majority of nearby 
ULIRGs are energetically dominated by compact (<500 
pc), highly dust-obscured nuclear cores, rather than by 
extended (>kpc), weakly obs cured (Ay < 15 ni ag) star- 
burst activity in host galaxies l|Soifer et alJ2000(l . There- 
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fore, a primary goal in the study of nearby ULIRGs is to 
establish whether the obscured compact cores are very 
compact starbursts or AGNs. Optical spectroscopy is 
an often used way to distinguish an AGN from a star- 
burst. In the presence of an AGN obscured by dust in 
a torus geometry, clouds that are photo-ionized by the 
radiation from the AGN develop along the torus axis 
above a torus s cale height (the s o-called narrow line re- 
gions: NLRs) l|Antonuccilll993|) . The NLRs produce 
optical emission lines whose flux ratios differ from the 
clouds in starburst galaxies; galaxies that contain an 
AGN obscured by torus-sha ped dust are classifi ed opti- 
cally as Seyfert 2s jVcillcux fc Osterbroci3ll987|) . Hence, 
this type of obscured AGN can be identified by using 
optical spectroscopy. However, the nuclear regions of 
ULIRGs contain sufficiently highly concentrated molec- 
ular gas and dust that putative AG Ns may be obscured 
by d ust in virtually all directions l|Sanders fc Mirabel 
I1996D . Such buried AGNs ^ without significant NLRs are 
classified opti cally as "non-Seyfert s" , and so very dif- 
ficult to find l)Maiolino et al J [20031) : however, they are 

^ We use the term "obscured AGN" if the AGN is obscured along 
our line of sight, and "buried AGN" if the AGN is obscured along 
virtually all sightlines. 
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predi cted to be abundant in the universe ijFabian et al.l 
l2002j) . In fact, the initial evidence of powerful 
buried AGNs has b een detected in several optical 
non-Seyfert ULIRCxS llP^idlev fc Wvnn-Williamsl [19971 
^ifcr et al. 2000: Imanishi & Dudlev '2000'; "Tran et al.' 
2001.: Imanishi et al 200 1: I manishi fc Malonev ,2003: 
Imanishi et al.ll2003l l200d[l . so that it is essential to un- 
derstand the contribution from buried AGNs to the total 
energy output of the ULIRG population. 

One of the most effective ways to distinguish a buried 
AGN from a very compact starbust as the primary en- 
ergy source in the core of a ULIRG is to investigate 
whether there is a strong X-ray-emitting source. A lu- 
minous buried AGN should exhibit intrinsically strong 
X-ray emission at E > 2 keV from the vicinity of the ac- 
cretion disk around the central supermassive black hole 
l)Elvis et al.lll994|) . whereas in a starburst, the intrin- 
sic X-ray emission above E' ~ 2 keV, rel ative to bolo- 
metri c luminosity, is significantly weaker ijRanalli et al.l 
12003"). Therefore, the detection of intrinsically strong, 
but highly absorbed. X-ray emission can provide strong 
evidence for a powerful buried AGN. Although mod- 
erately X-ray absorbed, Compton thin (Nh < 10^'* 
cm~^) AGNs are easy to find with currently available 
X-ray data at E = 2-10 keV fFranccschini e t alJl200l 
iPtak e t al. 2003; Imanishi & Tcrashima 2004), it is very 
difficult to unveil and quantify the luminosities of Comp- 
ton thick (Nh > lO^'' cm~^) buried AGNs using such 
data. Given that most of th e AGNs detected so far 
in optical Seyfcrt-2 ULIRGs fFranceschini et al.' '20Q^ 
[iwasawa et al, 2001; Braito ct al. 2003; Iwasawa et aLl 
120051) and non-Seyfert galaxies wi th Lir ^ 
llDone et all 119961: iVignati et a"niT999; Delia Ceca et all 
l2002^ suffer from Compton-thick absorption, it is very 
likely that the majority of the putative buried AGNs 
in optical non-Seyfert ULIRGs are also Compton thick. 
High-sensitivity observations at E > 10 keV using a next- 
generation X-ray satellite are required to properly eval- 
uate the energetic importance of buried AGNs in a large 
ULIRG sample. 

Another approach to investigating whether a ULIRG 
core has an intrinsically strong X-ray emitting source 
is to search for signatures of the effects of the lu- 
minous X-ray source on the surrounding interstellar 
medium. Around a buried AGN, X-ray dissociation re- 
gions (XDRs; Maloney et al. 1996) should develop, 
rather than the photo-dissociation regions (PDRs) usu- 
ally found in starburst regions. XDRs and PDRs can be 
distinguished if they show different flux ratios of emission 
lines. In particular, emission lines in the far-infrared to 
millimeter wavelength range are potentially a powerful 
tool, because low dust extinction at these wavelengths 
makes it possible to uncover the signatures of XDRs 
around powerful buried AGNs that are usu ally deeply 
embedded in dust. IKohno et all l|2001ll200^ and lKohnd 
(|20O5) proposed that the HCN (J=l-0) (A=3.38 mm), 
HCO+ (J=l-0) (A=3.36 mm), and CO(J=1~0) (A=2.60 
mm) lines are useful for this purpose, and argued that 
AGN-dominated galaxies (XDRs) show systematically 
larger HCN/HCO+ and HCN/CO ratios in brightness 
temperature than sta rburst-dominated galaxies (PDRs) . 
Ilmanishi et"an lj2004D pointed out that the HCN/HCO+ 
ratio is particularly effective for investigating the nature 
of heavily obscured ULIRG cores for the following rea- 



sons: (1) since both molecular lines probe high-density 
molecular gas (riHa >10'* cm^'^), their relative ratios are 
not strongly contaminated by more extended lower den- 
sity gas, whose relative fraction can vary among differ- 
ent galaxies, and (2) the wavelengths of the two lines 
are sufficiently close that they are observable simultane- 
ously with wide-band receiving systems, reducing the net 
observing time and systematic errors. Since HCN and 
HCO"'" are both moderately strong lines, this method 
is applicabl e to a l arger number of galaxies than that 
proposed bv lUsero~et al. ( 2004, ). which uses various faint 
molecular lines for the purpose of distinguishing between 
XDRs and PDRs. 

Interferometric maps are preferable to single-dish 
large-aperture observations for investigating the cores 
of ULIRGs in the millimeter wavelength range, be- 
cause they reduce the contamination from extended star- 
formation emission in the host galaxy. We applied this 
millimeter interferometric method to the infrared lumi- 
nous galaxy NGC 4418 (Lir, ~ 10"L(7>), wh ich likely con- 
tains a powerful bu ried AGN ijDudlev fc W ynn- William^ 
IT997tlSpoon et alj r2001). and were successful in finding 
sign atures of a buried A GN from the HCN/HCO+ ra- 
tio ijimanishi et aljr2004() . demonstrating the potential 
of this method. In this paper, we report the applica- 
tion of this interferometric HCN/HCO+ method to three 
ULIRGs: UGC 5101, Mrk 273, and IRAS 17208-0014. 
Throughout this paper, we adopt Hq = 75 km 
Mpc-\ nu = 0.3, and Ha = 0.7. 

2. TARGETS 

The target ULIRGs were selected based on the follow- 
ing criteria: 1) Strong HCN emission is detected with 
single-dish radio telescopes, so that the study of HCN 
and HCO+ emission-line ratios at the nuclei is feasible 
with the Nobeyama Millimeter Array (NMA), our ob- 
serving facility. (2) They are nearby (z < 0.06), so that 
redshifted HCN and HCO+ emission lines are observable 
with the NMA receiving systems. The three ULIRGs 
UGC 5101, Mrk 273, and IRAS 17208-0014 were se- 
lected as the first targets. Table 1 summarizes the de- 
tailed information about these sources. An angular scale 
of 1" corresponds to a physical size of 0.7-0.8 kpc for the 
three ULIRGs. 

UGC 5101 (z = 0.040) is included in the origi- 
nal ULIRG sample (Sanders_et al. 1988a) and is a 
well-stud ied ULIRG. I t is classified optically as a 
LINER l|Veilleux et al.l '1995D or as a Seyfert 2 in a 
slightly different c lassification s cheme by a d i fferent 
group ("Goncalves et al. 199^. 'Im anishi et"al] 1)200 If) 
and Imanishi & Maloney (2003) performed infrared 3- 
4 /im spectroscopy and detected polycyclic aromatic 
hydrocarbon (PAH) emission at 3.3 /xm, a starburst 
indicator, as well as signatures of a powerful AGN 
deeply buried at the core. They inferred that the in- 
frared luminosity of UGC 5101 is dominated by the 
buried AGN, rather than the dete cted starburst ac- 
tivity. Radio VLB I observations ijSmith et al.l 119981: 
iLonsdale et alJ 120031^ and a near-infr ared high-spatial- 
resolution image f Scoville et al J 120001) showed compact, 
high surface-brightness emission, which was interpreted 
as an AGN by these authors. Similar comp act emission 
was detected in the infrared at 12.5 /im ijSoifer et alJ 
Hmh . The presence of a powerful obscured AGN 
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was con firmed by subsequen t 2-10 keV X-ray ob ser- 
vations (|Imanishi et alJ 1200.^1 . lArmus et all l|2004D de- 
tected weak high-excitation forbidden emission Hnes orig- 
inating in the NLRs in a mid-infrared 5-40 /zm high- 
resolution spectrum. Judging from the d etection of NLR- 
originating emission hnes in the optical llGoncalves et al.l 
119991) and mid- infrared ijArmus et al.ll2004(l . there may 
be a small angular range that is transparent to the AGNs' 
ionizing photons (mostly UV and soft X-rays). HCN 
(J=l -0) da ta based on a single-dish telescope are avail- 
able ijGao fc Solomon ,,2004) and indicate that if a signif- 
icant fraction of the detected HCN emission comes from 
the nucleus of UGC 5101, then an investigation of the 
HCN/HCO+ ratio at the core should be feasible using 
the NMA with a reasonable integration time. 

Mrk 273 (z = 0.03 8) is also a well-studied ULIRG 
ijSanders et alJ Il988al). It is classified optically as a 
Sevfert 2 ijVeilleux et all 1 1999(1 . Its Seyfert-2 classifica- 
tion and the detection of strong NLR-originating high- 
excitation forbidden emission lines in the mid-infrared 
5-11 /im spectrum (Gcnzcl ct al. 199^ suggest that the 
dust distribution around the putative AGN is torus- 
shaped and that a significant amount of ionizing radi- 
ation is escaping the AGN. 2-10 keV X-ray data have re- 
vealed the presence of a moderately absorbed Compton- 
thin AGN (Iwasawa 1999; Xia et ai. 2002; Balestra ct al. 
I2005j) . The presence of an obscured AGN in Mrk 273 is 
also suggested from (1) a smaller equivalent width of 3.3 
yum PAH emission than the typic al value for starburst 
galaxies ijlmanishi fc Du dlev''200(7) . (2) deep 9.7 fj,m sili- 
cate dust absorption and weak PAH emission features in 
the 8-13 ^lul spectrum (Dudley 1999; Soifcr et al. 2002), 
and (3) compact hi gh surface brightn ess mid-infrared 
emission at 12.5 um ^Soifer et alJl200nD . The HCN fiux 
from Mrk 273 observed with a single-dish telescope is as 
strong as that from UGC 5101 ( Gao & Solomon 2004), so 
that investigations of nuclear HCN and HCO+ emission 
are also feasible for Mrk 273. 

IRAS 17208-00 14 (z = 0.042) is cl assified optically 
as an Hll-region ((Veilleux et all 1199,51) . Unlike other 
ULIRGs, the mid-infrared 8-23 ^m emission of IRAS 
1720 8—0014 is spatiall y extended, with no prominent 
core ijSoifer et al.ll200n() . Its radio emission as detected 
with the VLA is rel atively compact l)Martin et alJlT989t 
[Condon et al."1996^, but its surface brightness is not 
high enough to indicate strongly the presence of an AGN 
l|Diamond et al.lll99^ . Its 8-13 /im spectrum is typical 
of sta rburst galaxie s, with prominent PAH emission fea- 
tures l|Dudlevll 1999(1 Its 3-4 /xm spectrum also displays a 
clear 3.3 /im PAH emission feature, but the detection of a 
moderately strong 3.1 /im ice absorption feature suggests 
an energy source more centrally-concentr ated than the 
surrounding dust, su ch as a buried AGN ijRisaliti et al.l 
l2?M ITma,nis"hill2003) . However, in IRAS 17208-0014, 
the AGN signature is much weaker than UGC 5101 
and Mrk 273. The detected HCN fiux using a single- 
dish te lescope is the largest among the three ULIRGs 
lIcT^TTfc Solo mon 200l) . 

3. OBSERVATIONS AND DATA REDUCTION 

Interferometric observations of HCN (J=l-0) (Ajost — 
3.3848 mm or i^,.est ^ 88.632 GHz) and HCO+ (J=l-0) 
(Arest — 3.3637 mm or i^rest = 89.188 GHz) were per- 
formed with the NMA and RAINBOW Interferometer 



at Nobeyama Radio Observatory (NRO) between De- 
cember 2004 and March 2005. Table 2 summarizes the 
detailed observing log. The NMA consists of six 10-m 
antennas, and observations were performed using the D 
(the longest baseline was 82 m), C (163 m), and AB 
(351 m) configurations. The RAINBOW interferometer 
is a seven-element combined array that includes six 10- 
m antennas (NMA) and the NRO 45-m telescope. The 
RAINBOW observations were performed when the NMA 
array was in the AB configuration, and the longest base- 
hne was 410 m. The sensitivity of RAINBOW is about 
twice that of the NMA array in the 3-mm wavelength 
range, because the inclusion of the NRO 45-m telescope 
increases the total aperture size by a factor of four. 

The backen d was th e Ultra- Wide-Band Correlator 
(UWBC) ( Ok umura et a l. 2000), which was configured 
to cover 1024 MHz with 128 channels at 8-MHz reso- 
lution. The central frequency was set to 85.49, 85.67, 
and 85.24 GHz for UGC 5101, Mrk 273, and IRAS 
17208—0014, respectively, so that both the redshifted 
HCN and HCO^ lines were observable simultaneously. 
A bandwidth of 1024 MHz corresponds to ~3500 km s~^ 
at i/ ^ 86 GHz. The field of view at this frequency is 
~26" (full width at half maximum; FWHM) and -77" 
(FWHM) for RAINBOW and NMA, respectively. Since 
the Banning window function was applied to reduce side 
lobes in the spectra, the effective resolution was widened 
to 16 MHz or 55 km s"^ at ~ 86 GHz. To calibrate the 
passband across the 128 channels, the bright quasars 3C 
84 (UGC 5101), 0355+508 (Mrk 273), and 3C 454.3 and 
3C279 (IRAS 17208-0014) were observed every observ- 
ing day. Quasars 0954-f658, 1418-1-546, and 1741-038 
were used to calibrate temporal variation in the visibil- 
ity amplitude and phase for UGC 5101, Mrk 273, and 
IRAS 17208-0014, respectively. 

The UVPRO C-II package developed at NRO 
l|Tsutsumi et all Il99'7() and the AIPS package of 
the National Radio Astronomy Observatory were used 
for standard data reduction. Corrections for the antenna 
baselines, bandpass properties, and the time variation of 
the visibility amplitude and phase were applied to all of 
the data. A fraction of the data with large phase scatter 
due to bad radio seeing was removed from our analysis. 
After clipping a small fraction of data of unusually 
high amplitude, the data were Fourier-transformed 
using a natural UV weighting. The flux calibration was 
made using the observations of the appropriate quasars, 
1749-1-096 (UGC 5101), 1741-038 and 0355-1-508 (Mrk 
273), or 1741-038 (IRAS 17208-0014), whose flux 
levels relative to Uranus or Neptune had been measured 
at least every month between December 2004 and March 
2005. A conventional CLEAN method was applied to 
deconvolve the synthesized beam pattern. The total 
net on-source integration times were ^^29 hours (UGC 
5101), -8 hours (Mrk 273), and ~18 hours (IRAS 
17208-0014). The synthesized beam sizes are 5'.'8 x 
4'.'5 (position angle is 45° west of north; PA = —45°) for 
UGC 5101, 1'.'9 X 1'.'5 (PA = -66°) for Mrk 273, and 
3'.'1 X 2'.'4 (PA = -1°) for IRAS 17208-0014. Absolute 
positional uncertainties of NMA/RAINBOW maps were 
estimated to be much less than 1 arcsec. 

4. RESULTS 
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HCN (J=l-0) or HCO+ (J=l-0) data are available 
for several ULIR Gs (Solomon ct al. 199^ iCurran et al.l 
l2000t iGao fc SolomoDi i2004^) . but observations have 
generally been perfo rmed with single-dish telescopes. 
iRadford et al.' (1991*) showed an interferometric HCN 
map of the ultraluminous infrared galaxy Arp 220. With 
this exception, no detailed morphological information for 
HCN emission exists for ULIRGs. Therefore, this work 
allows us to reveal for the first time the spatial distri- 
bution of HCN (J=l-0) in the ULIRG population other 
than Arp 220. Ours are also the first published interfer- 
ometric maps of ULIRGs for HCO"'" (J=l-0) emission. 

Figure 1 shows the contours of the continuum emis- 
sion from UGC 5101, Mrk 273, and IRAS 17208-0014. 
Signals away from the HCN and HCO+ emission lines 
were combined to produce the continuum maps. Strong 
continuum emission is found in all ULIRGs and the con- 
tinuum peaks are spatially coincident with the nuclear 
positions of these ULIRGs. Table 3 summarizes the esti- 
mated continuum flux levels. The detection significance 
is > 6cr (UGC 5101), > 4cr (Mrk 273), and > IOct (IRAS 
17208—0014). The spatial extents of the continuum emis- 
sion agree with the beam sizes within <50%. 

After the continuum emission was subtracted, we in- 
vestigated the spatial distribution of the HCN and HCO"^ 
emission. Figure 2 shows the integrated intensity maps 
of the HCN emission for UGC 5101, Mrk 273, and IRAS 
17208-0014. Figure 3 displays the HCO+ integrated 
intensity maps for UGC 5101 and IRAS 17208-0014. 
HCO+ emission is not clearly detected from the nucleus 
of Mrk 273 in our interferometric map (< 3ct). With 
this exception, the HCN and HCO+ emission lines are 
clearly (> lOtr) detected. The detected emission lines 
are coincident with the nuclear positions of the ULIRGs 
and has spatial extents consistent with the beam sizes 
within <50%. For UGC 5101 and IRAS 17208-0014, 
the signal-to-noise ratios around the HCN and HCO+ 
lines are good enough to present channel maps, which 
are shown in Figures 4 and 5, respectively. 

Figure 6 shows continuum-subtracted nuclear spectra 
around the HCN and HCO^ emission lines. The pres- 
ence of HCO+ emission from Mrk 273 is unclear, but 
HCN emission from all the ULIRGs, and HCO+ emission 
from UGC 5101 and IRAS 17208-0014 are clearly recog- 
nizable. Figure 7 presents the Gaussian fits to these de- 
tected emission lines. The central velocity and line width 
of the Gaussian were determined independently for the 
HCN and HCO+ lines. The HCN emission from UGC 
5101 suggests the presence of double peaks, which may 
be related to the nuclear rotating molecular disk probed 
by the CO (J=l-0) line fGenzel et al. 1998). Therefore, 
a double Gaussian fit was also attempted for the HCN 
hue of UGC 5101. For the HCN and HCO+ hues from 
IRAS 17208—0014, we added a constant continuum to 
subtly adjust the residual continuum level in the vicinity 
of the HCN and HCO+ lines. Table 4 summarizes the 
Gaussian fitting results. 

The local standard of rest (LSR) peak velocities of 
HCN emission from Mrk 273 (11310 km s'^) and IRAS 
17208-0014 (12860 km s'^) agree, within 25 km s-\ 
with those of CO peak velocities (11324 km s""'^ for Mrk 
273 and 12837 km s"^ for IRAS 17208 -0014; Downes 
fc So lomon 1998). HCN and CO ( Downes fc SolomonI 
Il998j) line profiles are also similar for the nuclei of Mrk 



273 and IRAS 17208-0014. For the double Gaussian fit 
of HCN emission from UGC 5101, the velocity difference 
of the two components is ~400 km which is roughly 
comparable to that of CO (~480 km s~^; Genzel et al. 
1998). 

The integrated fluxes of HCN and HCO+ at the peak 
position, estimated from the Gaussian fits, are summa- 
rized in Table 5. For the HCN emission from UGC 5101, 
the single and double Gaussian fits give fluxes that are 
consistent to within <5%. Therefore, we adopt the val- 
ues based on single Gaussian fits for all sources. The 
derived HCN fluxes correspond to -60% (UGC 5101), 
-30% 2 (Mrk 273), and -55% (IRAS 17208-0014) of 
the HC N flux measurements f rom single-dish radio tele- 
scopes lIGao fc Solomonll200l . 

The HCN/HCO"'" ratios in brightness temperature (oc 
X flux density) are also summarized in Table 5. For 
the purpose of detecting the effects of putative AGNs 
on the surrounding interstellar medium at ULIRG cores, 
we regard our HCN/HCO^ ratios as a more appropri- 
ate tool than the ratios pre viously derived fro m single- 
dish telescope obser vations l)Nguven-0-Rieu et al. 1992^ 
iSolomon et alJlT99^ . because (1) we can reduce the con- 
tamination from the extended flux of the host galaxy, 
and (2) in our NMA data, the HCN and HCO+ fines 
were observed simultaneously with the same receiver and 
same correlator unit, so that possible systematic errors 
should be minimal. Although the absolute flux calibra- 
tion uncertainty of the NMA and RAINBOW data might 
be as large as —20%, this uncertainty does not affect 
the HCN/HCO+ ratios. The uncertainty in the ratios 
is dominated by statistical noise and fitting errors (see 
Figure 7). 

5. DISCUSSION 

5.1. Comparison of HCN/HCO^ Ratios with Other 
Galaxies 

Figure 8 shows the HCN/HCO+ and HCN/CO ra- 
tios in brightness temperature for UGC 5101, Mrk 273, 
and IRAS 17208-0014, together with previously ob- 
tained data for nearby starburst and Seyfert galaxies 
(Kohno ct al. 2001; Kohno 2005), and the AGN-powered 
infrar ed luminous galaxy NGC 4418 ijimanishi et alJ 
|2Q^). The ratios for the three main nuclei of the in- 
frared luminous merging galaxy Arp 299 (Imanishi et 
al. 2006, in preparation) are also plotted. iKohno et alJ 
f2001) and Kohno ( 2005J found empirically that AGN- 
dominated galaxies tend to have higher HCN/HCO"*" 
and HCN/CO ratios than starburst-dominated galax- 
ies: AGN-dominated galaxies should be distributed on 
the upper right-hand side, whereas starburst-dominated 
galaxies should occupy the lower left region in Figure 8. 

The HCN/CO ratio in brightness temperature (ab- 
scissa) may increase if the HCN abundance relative 
to CO is enhanced by X -ray radiation from an AGN 
(jLepp fc Dalgarnol Il996() . However, another possi- 

^ This value seems low, but we have no clear answer. The flux 
difference between the single-dish and interferometric observations 
may be due to discrepancies in the absolute flux scale between the 
two observations, or the presence of extended HCN emission unre- 
covered by our interferometric observations (less likely), or possi- 
bly a substantial overestimate of the HCN flux by Gao & Solomon 
(2004), coming from relatively large scatters in their single-dish 
data (their Figure 2). 
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ble explanation for the enhanced HCN/CO brightness- 
temperature ratio is an increase in molecular gas den- 
sity. Since the dipole moment of HCN ~ 3 debye; 
Millar et al. 1997) is much larger than that of CO 
(jj, ~ 0.1 debye; Millar et al. 1997), HCN emission 
traces a much higher density molecular gas (riu, > 10** 
cm~^) than is traced by CO (nna ^ 10^~^ cm"^). The 
HCN/CO brightness-temperature ratio can increase if a 
galaxy contains a larger fraction of high-density molecu- 
lar gas (e.g., Solomon et al. 1992; Heifer & Bhtz 1993; 
Kohno et al. 1999; Curran et al. 2000). In fact, ULIRGs 
are expected to have a larger fraction of high-density 
molecular gas than less infrared-luminous galaxies do 
fSandcrs & Mirabel 199^, so that HCN/CO brightness- 
temperature ratios can be high in ULIRGs, even without 
the HCN abundance enhancement caused by powerful 
AGNs. Furthermore, the HCN and CO data were taken 
on different observing runs, with different beam sizes, 
under different weather conditions, all of which could in- 
troduce additional systematic uncertainties. Therefore, 
the interpretation of HCN/CO brightness-temperature 
ratios in ULIRGs is not simple. 

Such uncertainties are much reduced if we use the 
HCN/HCO+ ratios. Both HCN and HCO+ {n ~4 debye; 
Botschwina et al. 1993) probe high-density molecular gas 
(riHa > lO'' cm"^), so that the HCN/HCO+ brightness- 
temperature ratio is insensitive to the different fraction 
of dense molecular gas, relative to diffuse one, in indi- 
vidual galaxies. Furthermore, since both the HCN and 
HCO^ emission lines were observed simultaneously with 
NMA/RAINBOW for all the sources plotted in Figure 
8, systematic errors are not a serious concern. For this 
reason, we use the HCN/HCO"'' ratios (ordinate) in our 
discussion of the possible presence of powerful AGNs. 

Our new millimeter interferometric data for UGC 5101 
and Mrk 273 reveal that the HCN/HCO+ brightness- 
temperature ratios are high and in the range expected 
for AGN-dominated galaxies (HCN/HCO+ > 2; Kohno 
2005). This is consistent with the fact that previous data 
have found the signatures of powerful AGNs in UGC 5101 
and Mrk 273 (§2). 

The interpretation of our data for IRAS 17208-0014 
is rather unclear. No strong evidence for a luminous 
buried AGN in IRAS 17208-0014 has been found at 
other wavelengths (§2). The HCN/HCO+ brightness- 
temperature ratio for IRAS 17208-0014 (--1.7) is located 
in between the expected values for AGN- and starburst- 
dominated galaxies. A powerful AGN possibly hidden in 
IRAS 17208-0014 may slightly increase the HCN/HCO+ 
ratio, but the presence of an AGN is not clear from our 
data. 

5.2. Interpretations of HCN/HCO+ Ratios 

In ULIRGs, the profiles of HCN (J=l-0) and CO 
(J=l-0) lines are generally similar (Solomon et al. 1992; 
Gao & Solomon 2004; This paper), despite the differ- 
ence of their optical depths. The similar line profiles are 
naturally explained b y the widely accepted mist model 
ijSolomon et al.llT98'^ . where a molecular gas cloud con- 
sists of a large number of small clumps, with a filling 
factor of <1 at any fixed velocity. In this model, even if 
each clump is optically thick for a molecular line, we can 
probe clumps at the background side of the molecular 
gas cloud. 



If each clump is optically thin for a molecular line, 
its flux will increase linearly with increasing abun- 
dance. However, each clump is thought to be opti- 
cally thick for CO (J=l -0) lines ijSolomon et alJll987t 
IScoville fc Sand ers' 1987^. It is not clear whether each 
clump is optically thick or thin for HCN (J=l-0) and 
HCO+ (J=l-0) lines, but we can roughly estimate their 
optical depths (r), based on the formula of 

TocA X X abundance^ (1) 

where A is the Einstein coefficient and is propor- 
tional to /i^, where n is the value of dipole moment 
(|Scovil le fc Sa nderslll987|) . If we adopt the HCN and 
HCO+ abund ance , relative to CO, are 10"** ~ 5x10"'^ 
ijBlake et al.l il987|) . the optical depths of both HCN 
(J=l-0) and HCO+ (J=l-0) lines are 0.2-2 times as 
large as that of CO (J=l-0) line. Thus, this simple es- 
timate suggests that each clump can be optically thick 
for HCN (J=l-0) and HCO+ (J=l-0) lines. Nguyen-Q- 
Rieu (1992) also showed that each clump can have optical 
depth slightly larger than unity for HCN (J=l-0) line, if 
the volume filling f actor of the clump is as small as f < 
0.1 HSolomon et al.lll987llKoda et al.ll2n(i5|) . 

For an optically thick case, although the relation be- 
tween flux and abundance is not so simple, their positive 
correlation is still expected for the following main rea- 
son. Assuming that each clump has a decreasing radial 
density profile (cx r~^-^; Gierence et al. 1992), line emis- 
sion comes ma inly from the layer w here an optical depth 
reaches unity (|Gierens et al.lll992|) . When a molecular 
abundance increases (decreases), the r = 1 layer for that 
molecular line moves outward (inward). Consequently, 
the surface area of the layer, or its area filling factor, be- 
comes larger (smaller) , increasing (decreasing) the flux of 
that molecular line. Unless the optical depth is extremely 
large, this effect can be significant and flux will increase 
with increasing abundance, if not completely linearly, 
as was demonstrated for the o ptically thick CO (J=l- 
0) line llVerter fc Hodgdll995t IWilso^lT995t ISakamotol 
'1996'; 'Arimot o et al.lll996D . ' Hence, qualitatively, the 
high HCN/HCO"'" ratio in the brightness temperature 
can be explained by an enhanced HCN/HCO+ abim- 
dance ratio in dense molecular gas from which the bulk 
of the observed HCN and HCO"'" fluxes originate. 

Theoretical models allow the calculation of the abun- 
dance of molecules, including HCN and HCO+, as a func- 
tion of distance from a central UV illuminating source 
or X-ray e mitting AGN deeply bu r ied in molecular gas 
and du st iLepp fc Dalgarnol Il996t iMeiierink fc Spaan^ 
12005ft . iLepp fc Dalgarnol l)199^ ound that the abun- 
dance of HCN can be enhanced relative to HCO+ 
if a strong X-ray source is pr esent inside molecular 
gas. IMeiierink fc Spaaii^ l|2005|) calculated the effects 
of a central X-ray-illuminating AGN on the surrounding 
molecular gas. They found that the HCN/HCO+ abun- 
dance ratio in uniformly-distributed dense molecular gas 
(lO^-^ cm~3)^ as probed by HCN and HCO+, around 
a luminous X-ray-emitting source can be substantially 
higher than the ratio in dense molecular gas illuminated 
only by a central UV source (their Fig. 10, model 4). 
Thus, we have some theoretical basis for the enhanced 
HCN abundance, relative to HCO+, in the presence of 
an X-ray-emitting AGN. 

However, further sophisticated model calculations 
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are definitely required to quantitatively interpret our 
observational results. Models available to date 
ijLePD fc Dalgarnolll 996l'M eiierink fc Spaansll2005|) have 
only calculated the HCN/HCO"'" abundance ratio and 
showed that the abundance ratio is a function of distance 
from the central illuminating source. To convert from the 
HCN/HCO"*" abundance ratio to brightness-temperature 
ratio, we have to know the region from where the bulk of 
the observed HCN and HCO+ fluxes come, which has not 
been explicitly mentioned in any available model calcula- 
tions. Furthermore, the models assume uniform dust/gas 
distribution for simplici ty. If the distribution is clumpy, 
which is more realistic (jSolomon et al.lll987j) . the abun- 
dance ratio can change substantially. 

In order for our millimeter interferometric 
HCN/HCO+ method to help to detect buried AGN 
signatures, high-density molecular gas must be concen- 
trated near the core region where putative AGNs reside 
and can affect the gas. If high-density molecular gas 
is distributed in between double nuclei, for example 
as in NGC 6240 (Nakanishi ct al. 2005;) , then our 
method may no longer work. Furthermore, s hocks 
mav enhance HCO"'" emission URawhngs et a1.l l200(l 
120041) . If shocks driven by outflows or superwind from 
starbursts coexist with powerful buried AGNs, and if 
there is a large amount of high-density molecular gas 
that is affected by the shocks, then the HCN/HCO+ 
brightness-temperature ratio may decrease compared 
to that of a pure buried AGN. Thus, some fraction of 
the galaxies containing powerful AGNs can show low 
HCN/HCO+ ratios and may be missed when searching 
based on the HCN/HCO+ ratio. 

Despite these remained ambiguities, this millimeter 
interferometric method can empirically be an effective 
means of revealing the presence of deeply buried AGNs 
at ULIRG cores, because of the negligible dust extinction 
in this wavelength range. If this method is established 
from a larger number of sources, future observations with 
high spatial resolution and high sensitivity using ALMA 
are potentially a powerful probe for optically and even 
2-10 keV X-ray-elusive buried AGNs in ULIRG cores. 

6. SUMMARY 

We performed millimeter interferometric observations 
of the HCN (J=l-0) and HCO+ (J=l-0) lines of the 
three ultraluminous infrared galaxies, UGC 5101, Mrk 
273, and IRAS 17208-0014, using NMA/RAINBOW. 
These are the first interferometric maps of a ULIRG sam- 
ple for which the HCN (J = 1-0) and HCO+ (J = 1-0) 
molecular lines were observed simultaneously. We ex- 
tracted the HCN to HCO"'" brightness-temperature ra- 
tios at the cores of these ULIRGs, and compared the 
ratios with those previously estimated for less-luminous, 
AGN-dominated and starburst-dominated galaxies. We 
then looked for possible signatures of the influence of 
AGNs, strong X-ray emitters, on the surrounding inter- 
stellar medium at the nuclei. We found that: 

1. Continuum and HCN emission were clearly de- 
tected in all ULIRGs. HCO+ emission was found 
in UGC 5101 and IRAS 17208-0014. The spatial 
distribution of the emission, when detected, is coin- 
cident with the nuclear positions of these ULIRGs. 



ratios in brightness temperature, as found in AGN- 
dominated galaxies. Previous observations had 
suggested luminous obscured AGNs in UGC 5101 
and Mrk 273, and so our millimeter interferomet- 
ric data may be revealing the signatures of the ef- 
fects of a luminous X-ray-emitting AGN on the sur- 
rounding interstellar medium at the nucleus. 

3. The HCN/HCO^" ratio for IRAS 17208-0014 was 
in between the regions occupied by AGN- and 
starburst-dominated galaxies. There was no clear 
evidence for a hidden AGN in IRAS 17208-0014 
based on previously obtained data at other wave- 
lengths. The presence of a luminous buried AGN 
in IRAS 17208—0014 remains uncertain based on 
currently available millimeter interferometric data. 

4. Although the sample size is still small, this mil- 
limeter interferometric HCN/HCO+ method suc- 
ceeded in confirming AGN signatures in the two 
ULIRGs (UGC 5101 and Mrk 273) that were pre- 
viously known to possess luminous obscured AGNs. 
Since dust extinction is negligible at this wave- 
length range, we argue that our method is po- 
tentially an important tool for unveiling powerful, 
but deeply buried, AGNs in optical "non-Seyfert" 
ULIRGs in the ALMA era. 
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TABLE 1 

Detailed information about the observed ULIRGs 



Object 


Redshift 


/l2 


/25 


/eo 


/lOO 


log iiR (log Lir/Lq) 


Far-infrared 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(ergs s"-"-) 


Color 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


UGC 5101 


0.040 


0.25 


1.03 


11.54 


20.23 


45.5 (12.0) 


cool 


Mrk 273 


0.038 


0.24 


2.28 


21.74 


21.38 


45.7 (12.1) 


cool 


IRAS 17208-0014 


0.042 


0.20 


1.66 


31.14 


34.90 


45.9 (12.3) 


cool 



Note. — Column (1): Object. Column (2): Redshift. Columns (3)-(6): fi2, £25, feo, and fioo are the IRAS FSC Ruxes 
at 12, 25, 60, and 100 fim, respectively. Column (7): Decimal logarithm of the infrared (8—1000 fim) luminosity in ergs s~^ 
calcu lated as follows: Lir = 2.1 x lO^'^x D(Mpc)^ x (13.48 x /12 + 5.16 x /25 + 2.58 x /eo + /100) ergs s"^ JSanders fc Mirabell 
Il99tf) . The values in parentheses are the decimal logarithms of th e infrared luminositie s in units of solar luminosities. Column 
(8): ULIRGs with fi2/f25 < (>) 0.2 are classified as cool (warm) JSanders et al.lll988bf) . 



TABLE 2 
Observing log 



Object 




Array 
Configuration 




Observing Date 
(UT) 


UGC 5101 
Mrk 273 

IRAS 17208-0014 




D 

RAINBOW 
C 

RAINBOW 
AB 
C 




2004 Dec 22, 23, 24 

2005 Jan 20 
2005 Mar 9, 10 
2005 Jan 22 and 31 
2005 Jan 26-28, Feb 7 
2005 Mar 6, 7, 16-19, 26 


TABLE 3 
Continuum emission 


Object 




Flux 
(mjy) 




Frequency 
(GHz) 


UGC 5101 
Mrk 273 
IRAS 17208-0014 




6 
6 

9 




85.5 
85.6 
85.3 




TABLE 4 

Gaussian fitting parameters of HCN and HCO+ 


emission lines 




Object 
(1) 


HCN 

(2) 


LSR velocity 
(km s^-"-) 

HCO+ 

(3) 


HCN 

(4) 


FWHM 

(km s"-"-) 

HCO+ 
(5) 


UGC 5101 

Mrk 273 
IRAS 17208-0014 


11830 11990 
11585+11985 11990 
11310 

12860 12820 


700 
270+375 
450 
530 


330 
330 

470 


Note. — Col.(l): 
velocity of the HCO+ 


Object name. Col. (2): 
emission peak. Col. (4): 


LSR velocity {vopt = (^ - 1) x 
Line width of the HCN emission 


c} of the HCN emission peak. Col. (3): LSR 
at FWHM. Col. (5): Line width of the HCO+ 



emission at FWHM. 
^ A double Gaussian fit of the HCN emission line (see §4). 
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TABLE 5 

HCN, HCO+, AND CO EMISSION 



Nucleus 


HCN 


HCO+ 


HCN/HCO+ 


CO 


Reference 




(Jy kms-l) 


(Jy km s~^) 




(Jy kms-l) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


UGC 5101 


6.6 


2.7 


2.5 


50 


a 




3.6+2.7 




>2,~1 






Mrk 273 


4.8 


<2.7 


>1.8 


78 


b 


IRAS 17208-0014 


14 


8.1 


1.7 


132 


b 



Note. — Col.(l): Object name. Col. (2): Integrated HCN intensity at tlie nucleus. Col. (3): Integrated HCO^ intensity at 
the nucleus. Col. (4): HCN/HCO"'' ratio in brightness temperature (cx x flux density) at the nucleus. The ratio is not affected 
by possible uncertainties in the absolute flux calibration (see §4). Col. (5): Integr ated CO intensity at the nu cleus taken from the 
literature. Col. (6); References for interferometric CO flux measurements, (a): IGenzel et alJ il998f) .: (b): IDownes fc SolomonI 
(|l99i). 



^ Blue (Left) and red (Right) components for the double Gaussian fit (see §4). HCO^ emission has a peak velocity close to 
that of the red HCN component. We attribute most of the HCO^ flux to the red component. 
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Fig. 1.— Continuum maps of UGC 5101, Mrk 273, and IRAS 17208-0014. The crosses show the coordinates of the ULIRG nuclei. 
The coordinates in J2000 are {9''35'"51.61= , +61°21'11'.'4) for UGC 5101, (13''44'"42.12^ +55°53'13('5) for Mrk 273, and (17''23™21.95=, 
-00°17'00'.'7) for IRAS 17208-0014. The nuclear coordinates of UGC 5101 are adopted from VLA radio data (Condon et al. 1991), while 
those of Mrk 273 and IRAS 17208—0014 are based on the CO peaks in millimeter interferometric maps llDownes &: Solomon 1998i) . The 
contours are —3, 3, 4.5, and 6 mjy beam~^ for UGC 5101, —4.5, —3, 3, 4.5, and 6 mjy boam^^ for Mrk 273, and —4, 4, 6, and 8 mjy 
beam-i for IRAS 17208-0014. 
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Fig. 2. Integrated intensity map of HCN emission. The contours arc 2.4, 3.6, 4.8, and 6.0 Jy km s ^ beam ^ for UGC 5101, 2.4, 3.2, 
and 4.0 Jy km s^^ beam^i for Mrk 273, and 4, 6, 8, and 10 Jy km s^^ beam^i for IRAS 17208-0014. 
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Fig. 3. — Integrated intensity map of HCO+ emission. The contours are 1.8, 2.7, and 3.6 Jy km s"! beam"! for UGC 5101, and 2, 4, 
and 6 Jy km s"^ beam"! for IRAS 17208-0014. 
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Fig. 4. — Channel maps of HCN and HCO+ emission for UGC 5101. Left: HCN emission. The contours arc —6, 6, and 10 mjy bcam"^. 
The r.m.s. noise level is ~2.5 mJy beam~^. Right: HCO+ emission. The contours are —6, 6, and 10 mJy beam~^. The r.m.s. noise level 
is ~2.5 mJy beam"^. 
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Fig. 5.— Charmcl maps of HCN and HCO+ emission for IRAS 17208-0014. Left: HCN emission. The contours arc ~9, 9, 15, 21, and 
27 mJy beam~^. The r.m.s. noise level is ~3 mJy beam^^. Right: HCO+ emission. The contours are —9, 9, 12, and 15 mJy beam~^. 
The r.m.s. noise level is ~3 mJy beam~^. 
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Fig. 6. — HCN and HCO^ spectra at the nuclei of the ULIRGs. The abscissa is the observed frequency in GHz and the ordinate is the 
flux in mjy beam~^. 
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Fig. 7. — Gaijssian fits to the HCN and HCO+ emission lines. The abscissa is the LSR velocity {vopt = — 1) x c} in km s ^ and 
the ordinate is the flux in mjy beam~^. Single Gaussian fits are used as defaults and are shown as dashed lines. For the HCN line of UGC 
5101, since there is evidence for a double peak, a double Gaussian fit is also attempted and shown as the dotted line. 



HCN/HCO+ in ULIRGs 



15 



o 
o 

+ 

O 

o 



ro - 



c\l - 



1 1 



pure AGN 



I I 




U5101 



■M273 
. O 



- (D)qOO 



N4418 
IR17208 








0.2 0.4 
HCN/CO ratio 



0.6 



Fig. 8.— HCN/HCO+ (or dinate) and HCN/CO (abscissa) ratios in brightness temperature (oc x flux). UGC 5101, Mrk 273, IRAS 
17208-0014, and NGC 4418 llmanishi et alj|2051) are plotted as filled stars. The three main nuclei of Arp 299 (Imanishi et al. 2006, in 
preparation) are also plotted as filled stars on the lower left side. The HCN/HCO"'' ratios are derived directly from our interferometric 
data, while the HCN/CO ratios are derived by combining our data with CO data in the literature. The plot of NGC4418 has been updated 
from Imanishi et al. 12004) to include new interferometric CO data by Dale et al.(2005). Other data points are taken from lKohno et all 
OOOU and .Kohna L2005.) . where sources with AGN-like (starburst-like) ratios are marked with filled squares (open circles). 



